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We report an experimental study of the properties of the velocity boundary layer in 
turbulent Rayleigh-Benard convection in a cylindrical cell. The measurements were made 
at Rayleigh numbers Ra in the range 2.8 x 10^ < Ra < 5.6 x 10^ and were conducted 
with the convection cell tilted with an angle 9 relative to gravity, at 6 — 0.5°, 1.0°, 2.0°, 
and 3.4°, respectively. The fluid was water with Prandtl number Pr = 5.3. 

It is found that at small tilt angles {9 ^ 1°), the measured viscous boundary layer 
thickness Sy scales with the Reynolds number Re with an exponent close to that for a 
Prandtl-Blasius laminar boundary layer, i.e. Sy ^ i?e~° ''^^° °^. For larger tilt angles, 
the scaling exponent of 6^ with Re decreases with 9. The normalized mean horizontal 
velocity profiles measured at the same tilt angle but with different Ra are found to have 
an invariant shape. But for different tilt angles, the shape of the normalized profiles is 
different. 

It is also found that the Reynolds number Re based on the maximum mean horizontal 
velocity scales with Ra as Re ~ iJa" "'^ and the Reynolds number Re^ based on the 
maximum rms velocity scales with Ra as Rea ~ Ra'^'^^, with both exponents do not 
seem to depend on the tilt angle 9. 

Several wall quantities are also measured directly and their dependency on Re are 
found to agree well with those predicted for a classical laminar boundary layer. These 
are the wall shear stress t (~ i?e^'^^), the viscous sublayer 6^ (~ i?e°-'^^), the friction 
velocity Ut i?e~°'^^) and the skin friction coefficient Cf i?e~°'^^). Again, all these 
near-wall quantities do not seem to depend on the tilt angle. 

We also examined the dynamical scaling method proposed bys Zhou and Xia [Phys. 
Rev. Lett. 104, 104301 (2010)] and found that in both the laboratory and the dynamical 
frames the mean velocity profiles show deviations from the theoretical Prandtl-Blasius 
profile, with the deviations increase with Ra. But profiles obtained from dynamical scaling 
in general have better agreement with the theoretical profile. It is also found that the 
effectiveness of this method appears to be independent of Ra. 



1. Introduction 

1.1. Rayleigh-Benard convection 

Rayleigh-Benard (RB) convection, which is a fluid layer heated from below and cooled 
from the top, is an idealized model to study turbulent flows involving heat transport and 
has attracted much attention during the past few decades ( Siggia|[T994 Kadanoff||2001 



Ahlers, Grossmann & Lohse 2009 Lohse & Xia 2010). The system is characterized by 



two control parameters; the Rayleigh number Ra, Prandtl number Pr, which are deflned 
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as 



and 



Ra 



igATH^ 



Pr = 



(1.1) 



(1.2) 



respectively. Here a is the thermal expansion coefficient, g the gravitational acceleration, 
AT the temperature difference between the bottom and the top plates, H the height of 
the fluid layer between the plates, v the kinematic viscosity, and k the thermal diffusivity 
of the converting fluid. In addition, the aspect ratio F = D/H [D is the lateral dimension 
of the system) also plays an important role in the structures and dynamics of the flow. 

In a fully developed Rayleigh-Benard turbulent flow, most of the imposed temperature 
difference is localized in the thermal boundary layers near the surface of the top and 
bottom plates, within which heat is transported via conduction ( Wu fc Libcharber|1991 



Belmonte, Tilgner k, Libchaber||1994 Lui fc Xia||1998 ). The velocity field has the same 



character: velocity gradient is localized in a thin layer near the plates, which is called vis- 
cous boundary layer. Turbulent flow in the central region of the RB cell is approximately 
homogenous and isotropic ( Zhou, Sun fc Xia|2008 Ni, Huang fc Xia|2011a 2012). As the 
top and bottom boundary layers contribute the main resistance to heat transfer through 
the cell and thus dominantly determine the Nusselt number, they deserve special atten- 
tion. Indeed, nearly all theories in RB convection are in essence boundary-layer (BL) 
theories. For example, a turbulent BL was assumed in the early marginal stability theory 



(Malkus 1954) and also in the models by Shraiman fc Siggia (1990) and [Siggia (1994) 
and by Dubrulle (2001 2002). On the other hand, a Prandtl-Blasius (PB) type laminar 



BL was assumption in the Grossmann & Lohse (GL) theory (Grossmann fc Lohse||2000 



2001 2002 2004). Therefore, direct characterization of the BL properties is essential for 



testing and differentiating the various theoretical models, and will also provide insight 
into the physical nature of turbulent heat transfer. 

1.2. Boundary layer measurements in turbulent thermal convection 

One of the earlier measurements of temperature and also velocity profiles in turbulent RB 
convection was taken by |Tilgner, Belmonte fc Libchaber ( 1993) in water {Pr = 6.6) at the 
fixed Ra = 1.1 x 10^ and at a fixed lateral position. Belmonte, Tilgner fc Libchaber ( 1993 ) 



extended these measurements over the range 5 x 10*^ ^ Ra ^ 10^^ in compressed gas 



(air) at room temperature (Pr=0.7), but still at fixed lateral position. Lui fc Xia (19981 



measured the mean temperature profiles at various horizontal positions on the lower 
plate of a cylindrical convection cell, the result shows that the thermal layer thickness 
5th varied over the plate for the same i?a, and the thinnest BL is closed to the center 



of the plate. Wang & Xia (2003) found similar results for a cubic cell, du Puits et al. 



(20076) measured high-resolution temperature profiles in Rayleigh-Benard convection 



near the top plate of a cylindrical container with air (Pr — 0.7) as the working fluid. 
Their result shows that the thermal BL thickness Sth ^ i?a~°'^''^ in the cell with F = 1.13. 



Sun, Cheung fc Xia ( 2008 ) found that the thermal boundary layer thickness scales with 
Ra~^ '^^ in a rectangular cell, at Pr = 4.3 and Ra ranging from 10* to 10^°. 

For the velocity measurement, the methods for determining the velocity profiles near 
the solid walls of the cell are developed in recent years. Since strong temperature fluc- 
tuations exist in Rayleigh-Benard convection, the well-established hot-wire anemometry 
could not be applied to this system. For the viscous boundary layer, the large temper- 
ature fluctuations make conventional laser Doppler velocimetry ineffective because the 
temperature fluctuations cause fluctuations in the refractive index of the fluid that in turn 
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make it difficult to steadily focus two laser beams to cross each other in the fluid (jXia, 



Xin fc Tong|p95| .[Tii gner et aL\ (11993) introduced an electrochemical labeling method 



and measured the velocity profile and boundary layer thickness near the top plate of a 



cubic cell filled with water, but only at a single value of Ra. In a later study, Belmonte, 



Tilgner & Libchaber (1993 1994) developed an indirect method — the correspondence 



between the peak position of the cutoff frequency profile of the temperature power spec- 
trum and the peak position of the velocity — to infer the viscous boundary boundary 
layer thickness in gaseous convection. This method has subsequently been used to infer 
the viscous layer in thermal convection in mercury (Naert, Segawa fc Sano|1997] ). It may 
just be that the method works in certain situations, but there is no theoretical basis for 
it. Xin, Xia & Tong (1996), using a novel light-scattering technique developed by Xia 



et al. ( 1995 ), conducted the first direct systematic measurement of velocity profiles in RB 
convection in a cylindrical cell as a function of Ra. They found S^, Ra~^'^^ from the 
velocity profile above the center of the lower plate. Qiu & Xia ( 1998a|& ) extended these 
measurements to convection in cubic cells, finding the same scaling exponent —0.16 at 
the bottom plate , but a t the sidewall a different result S^^ ^ i?a~°'^^. Using various or- 
ganic liquids. Lam et al. (2002j) explored the Pr dependence, finding S^, ^ Pr'^-^'^Ra~^-^^ . 
With the measured Ra — Re scaling relationship Re ^ Ra'^-^ obtained in these studies 
(in fact, it was the Peclet number Pe = vL/k rather than Re in some of these studies; 
please see Sun & Xia (2005) for more detailed discussions), the above results imply a 
scaling relation ~ /jg-u.ai: jj^ j-g^ent years, the technique of particle image velocimetry 



(PIV) has been introduced to the experimental study of thermal convection ( Xia, Sun 
fc Zhou|[2003l |Sun, Xi fc Xial|2005a[ |Sun, Xia fc "^^|20056[ |Sun fc Xia|[2()05[ ). |Sun 



et aL| ( |2008 ) further applied the PIV technique to study the viscous BL in a rectangular 
cell. Their results show that S^, ~ Ra~°'^'^ and S^, ~ Re~^'^, which showed that the 
viscous BL in thermal turbulence has the same i?e-scaling as a Prandtl-Blasius laminar 
BL. This result validates the laminar BL assumption made in the GL model in a scaling 
sense. Thus it appears there is a discrepancy in the measured scaling exponent of S^j 
with respect to Ra (Re) between those obtained in cylindrical and cubic cells and that 
obtained in rectangular cells. In Sun et al. (2008) it was argued that because of the more 
complicated fiow dynamics of the large-scale circulation (LSC), such as the azimuthal 
motion in the cylindrical cell (Sun et al. 20056 Brown et al. 2005 Xi et al. 2006) and 
the secondary flows in the cubic cell ( Qiu & Xia 1998a[ ), the shear flow near the plates 
is less steady as compared to that in the rectangular cell which is more close to quasi- 
two-dimension (quasi-2D). As the viscous boundary is created by the shear of the LSC, 
this may plausibly change the BL properties, resulting in a different exponent. However, 
the above argument has not been substantiated experimentally. Part of the motivation 
of the present work is to determine how the three-dimensional LSC dynamics will affect 
the BL properties. It is known that titling the cell by a small angle will "lock" the LSC 
in a fixed azimuthal plane in the sense that it will limit the range of the LSC's azimuthal 
meandering (Sun et al. 2005a) and reduce its azimuthal oscillation amplitude near the 



top and bottom plates of the cell (Ahlers, Brown & Nikolaenko 2006). In the present 
work, we present measurements of BL properties in a cylindrical cell with the cell titled 
over a range of angles. For small titling angles, we measure a boundary layer under a 
more steady shear comparing to the "leveled" case when the LSC can freely meander 
in the azimuthal direction but presumably the BL is otherwise unperturbed under such 
a small titling angle 1°). We also examine how the BL scaling exponent and other 
BL properties behave when the titling angle becomes not so small 1°), which would 
amount to a perturbation to the BLs. Boundary layers play such an important role in 
turbulent thermal convection, it is therefore important to examine how BLs respond to 
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external perturbations. Understanding the stability or instability of BLs is also relevant 
to the search for the so-called ultimate state of thermal convection, as the transition from 
the "classical state" to the ultimate one is essentially an instability transition of the BL 
from being laminar to being turbulent. 

In addition to the scaling of the BL thickness, the shape of the velocity profile near 
the top and bottom plates has attracted a lot of attention recently. Although the BL has 
been found scaling wise to be of Prandtl-Blasius type (at least in the quasi-2D case) , the 
time-averaged velocity profiles are found to differ with the theoretically predicted one 



(du Puits et al. 2007a Sun et al. 2008), especially for the region around the thermal 



BL. Recently, Zhou & Xia (20101 have proposed a dynamic scaling method that shows 



that the mean velocity profile measured in the laboratory frame can be brought into 
coincidence with the theoretical Prandtl-Blasius laminar BL profile, if it is resampled 
relative to the time-dependent frame that fluctuates with the instantaneous BL thickness. 
This method was tested initially for the case of velocity profile in turbulent convection 
in a quasi-2D rectangular cell with water as working fluid [Pr — 4.3). In a follow-up 



study using two-dimensional DNS data, Zhou et al. (2010) found that the method is also 
valid for thermal boundary layers and for the case of Pr = 0.7 as well. More recently, 
these authors further shown, again using numerical data, that the method works also 
in other positions in the horizontal plate other than the central axis dZhou et aL][20TT|) 
and in three-dimension (3D) cylindrical cell for moderate values of Ra { Stevens et al. 
20121). However, iScheel, Kim & Whitel (|2012|) and iShi, Emran & Schumacherl (|2012|), 



both using numerical approaches, have found that dynamic scaling works less well in the 
3D cylindrical geometry than in the quasi-2D case. However, the method has not been 
tested experimentally so far in a 3D system. Here we would like to examine the dynamical 
scaling method using the experimentally obtained instantaneous velocity proflles in our 
three-dimensional cylindrical cell. 

1.3. Organization of the paper 

The remainder of this paper is organized as follows. We give detailed descriptions of 
the experimental setup and measurement instrumentation in f|2] and present and analyze 
experimental results in f|3j which are divided into six subsections. In §3.1| we present the 
measured temperature profiles and corresponding position-dependent fluid properties, 
which will be used to calculate the viscous and Reynolds stresses. In j |3.2[ the measured 
velocity profiles and their characterizations are presented. In |3.3[ the scaling properties, 
with both Ra and Re, of the thickness dy obtained from the mean velocity profiles and 
6a obtained form r.m.s. velocity profiles, are presented and discussed. We also discuss 
the influence of the cell tilting angle 9 on the boundary layer scaling. In |3.4| statistical 
properties (r.m.s. and skewness) of the velocity field in the boundary layer region are 
discussed. In |3.5| we present results of the viscous and Reynolds shear stresses distri- 
butions in the boundary layer, and discuss the scaling of the wall quantities. In §3.6| we 
test the dynamic scaling method with respect to the measured instantaneous velocity 
profiles. We summarize our findings and conclude in ^ 



2. Experimental apparatus 

2.1. Convection cell 

The measurements were made in a cylindrical Rayleigh-Benard convection cell, which has 
been described in detail previously ( |Zhou fc Xia|2002[ [Sun et aL|2005&||Ni et aL|20116[ ). 
Here we give only its essential features. The top and bottom conducting plates are made 
of pure copper with a thin layer of nickel to avoid oxidation. The sidewall is made of 




Plexiglas. To avoid distortions in the images viewed by the camera, a square-shaped jacket 
is fitted around the sidewall of the convection cell. As shown in figure [T|b), the jacket is 
filled with water. The diameter and height of the cell is D = 19.6 cm and H — 18.6 cm, 
respectively. The aspect ratio F = D/H is thus close to 1. Two (three) thermistors are 
embedded in the top (bottom) plate. The top plate temperature is maintained constant 
by a refrigerated circulator (Polyscience Model 9702) that has a temperature stability 
of 0.01°C. A NiChrome wire (26 Gauge, Aerocon Systems) surrounded by fiberglass 
sleeving and Teflon tape is distributed inside the grooves carved under the bottom plate. 
The wire is connected with five DC power supplies (GE Model GPS-3030) in series to 
provide constant and uniform heating. During the measurement, the whole cell is placed 
in a homemade thermostat box that is kept at the same temperature (30°C) as that of 
the fluid at the centre of the cell. During the experiment the cell was tilted by an angle 6 
such that the circulation plane of the LSC was parallel to the image plane of the camera 
(the x-z plane, see flgure 1). 



2.2. PIV measurement 
The application of PIV to thermal turbulence has been described in detail in several 



previous publications (Xia et al. 2003 Sun et al. 2005 & 2008). Here we only provide 



details concerning the particular features of the present experiment. The PIV system 
consists of one GGD camera with 2048 x 2048 pixels, a dual pulse Nd-YAG laser with 
135 mJ per pulse, a synchronizer and software. As the cell was titled, both the GGD and 
the laser light-sheet were titled accordingly with the same angle. A 105 mm focal-length 
macro lens was attached to the GGD to achieve a measuring area with size varying from 
18 X 18 mm^ to 30 x 30 mm^. Each 2D velocity vector is calculated from a subwindow 
(32 pixels X 32 pixels) that has 50% overlap with its neighboring subwindows, so each 
vector corresponds to a region of 16 pixels x 16 pixels and each velocity map contains 
127 X 127 velocity vectors in the x-z plane (see figure 1). This corresponds to spatial 
resolutions of about 0.135 x 0.135 mm^ to 0.236 x 0.236 mm^ for velocities u and w 
measured in the horizontal x and vertical z directions, respectively. For the measurement 
&i 9 = 3.4° particles with diameter 2 /im were used, while particles with diameter of 
10 /im were used for measurements with other three tilted angles. For each run, typically 
about 25200 image pairs were acquired with frame rate of 2 Hz. 
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Figure 2. (a) A profile of mean temperature (T) measured at Ra = 6.8 x 10* and = 3.4°. 
(b) A profile of dynamic viscosity obtained from the mean temperature profile in (a). 



e 


Ra 


Pr 


Umax 


5v 


e 


Ra 


Pr 


Umax 


5v 


ieg.) 






(mm/ s) 


[mm) 


[deg.) 






[mm/ a) 


{mm 


0.5 


5.77 X 10* 


5.45 


4.77 


2.80 


2.0 


1.34 X 10" 


5.42 


7.47 


2.57 




2.79 X 10* 


5.41 


3.57 


3.40 




5.85 X 10* 


5.42 


5.21 


3.25 




1.55 X 10*^ 


5.39 


7.38 


2.44 




5.53 X 10" 


5.32 


13.49 


1.97 




2.93 X 10^ 


5.39 


9.11 


2.13 


3.4 


2.74 X 10* 


5.67 


4.02 


4.54 




4.26 X 10^ 


5.38 


10.88 


2.06 




4.19 X 10* 


5.66 


4.70 


4.23 


1.0 


1.68 X 10^ 


5.45 


7.53 


3.46 




6.69 X 10* 


5.66 


6.00 


3.58 




3.19 X 10^ 


5.39 


9.96 


2.87 




1.38 X 10" 


5.68 


7.59 


2.70 




5.54 X 10" 


5.42 


11.60 


2.66 




2.19 X 10" 


5.64 


9.40 


2.50 




9.46 X 10* 


5.45 


5.72 


3.82 




2.66 X 10" 


5.66 


10.20 


2.11 




2.40 X 10* 


5.42 


3.06 


5.07 




9.89 X 10* 


5.66 


6.47 


3.00 




6.00 X 10* 


5.43 


4.71 


3.94 




3.20 X 10" 


5.65 


10.84 


2.14 


2.0 


2.76 X 10* 


5.41 


3.76 


4.12 




4.28 X 10" 


5.63 


12.65 


1.89 




2.78 X 10" 


5.42 


10.09 


2.15 




5.19 X 10" 


5.56 


13.49 


1.84 



Table 1. Control parameters of the experiment: cell tilt angle 9, the Rayleigh number Ra 
and the Prandtl number Pr; and the measured maximum horizontal velocity Umax and viscous 
boundary layer thickness &v The data are listed in chronological order. 



3. Results and discussion 

PIV measurements were made at four values of the titling angle = 0.5°, 1°, 2°, 
and 3.4°. For each 9, measurements over a range of Ra were made. Table [T] lists the 
parameters {9, Ra and Pr) of each measurement, which typically lasted for about 3.5 
hours. As already mentioned, titling the cell by a small angle has the effect of "locking" 
the LSC's circulation plane at a fixed azimuthal angle (in reality it restricts the angular 
range of the LSC's azimuthal meandering). Thus, measurements made with small 9 are 
aimed at studying BL properties under more steady shear, but the BL itself is assumed to 
be unperturbed otherwise. For large values of 9 we wish to examine how the BL responds 
to relatively large perturbations. 

3.1. Temperature profile and fluid properties 

The local values of fluid properties are needed in calculating the viscous and Reynolds 
shear stresses, which requires measurement of the local temperature. Temperature profiles 



for the leveled case have been measured systematically by Lui & Xia ( 1998 ) in a similar 
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Figure 3. Coarse-grained vector maps of the instantaneous (a) and time-averaged (b) velocity 
field measured near the center of the bottom plate [Ra = 4.2 x 10* with 6 — 3.4°), the velocity 
scale bar is in unit of mm/s. 



cylindrical cell. To check whether titling the cell by a relatively large angle will change 
the temperature profile, we measured one mean temperature along the central axis {x = 
y = 0) of the cell at a titling angle 9 = 3.4° (Ra = 6.8 x 10^). The result is shown in 
figure ^a,) and the dynamic viscosity corresponding to the local temperature is shown 
in figure|2|b). As these results are similar to those obtained in previous studies by Lui 
k Xia| ( [1998^ and |Sun et al] ( |2008| , we will use results from those studies at similar Ra 
in the calculations of Reynolds stress (Sec. 3.5 ) and other wall quantities that require 
position dependent viscosity (density). 

3.2. Velocity profiles and the Reynolds number scaling 

Figure [3] (a) shows an example of measured instantaneous velocity map and (b) time- 
averaged velocity field taken over a period of 3.5 h (corresponding to 25200 velocity 
frames), with the cell tilted at9 = 3.4° and at Ra = 4.2x 10*. In the present measurement, 
X spans from —8.75 mm to 8.75 mm, and z spans from to 17.5 mm. From the velocity 
scale in figure [3] (a) and (b), it is seen that there exist velocity bursts with values much 
larger than the maximum velocity in the time-averaged velocity field. It is found that 
velocity maps measured at other tilt angles have similar features. As the mean velocity 
and the velocity fiuctuations do not exhibit any obvious dependence on the horizontal 
position X over the small range of the measurement, the quantities presented below are 
based on values averaged along the x-direction over the width of the measuring area. 

Figure |4] plots the velocity profiles for different tilt angles and various values of Ra, 
which shows that the shapes of the profiles are rather similar at this level of detail. 
Figurejsjplots normalized profiles in which U{z) is normalized by the maximum horizontal 
velocity Umax (Ra) (for ease of reference the values of Umax are also listed in Table [T]) 
and the distance z from the wall by the viscous boundary layer thickness Sy{Ra) (to be 
defined below) . The figure shows that up to 2^^, profiles for different Ra and for the same 
tilt angle collapse on to a single curve quite well (except perhaps those correspond to 
the largest Ra for 9 — 2.0 and 3.4°). Note that z ~ 2J„ is around where U reaches its 
maximum value and beyond this position it decays toward cell centre. So this position 
may be taken as the separation between the boundary layer region and the bulk. The 
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Figure 4. Time-averaged horizontal velocity profiles measured at tilt angles 9 — 0.5° (a), 1.0° 
(b), 2.0° (c) and 3.4° (d). In each plot the corresponding value of Ra decreases from top to 
bottom (see Table IT] for exact values). 



above results suggest that for the same tilt angle the profiles in the boundary layer region 
have an invariant shape with respect to different values of Ra. This result is consistent 
with the finding by Sun et al. (2008). In figure [s] we also plot the theoretical Prandtl- 
Blasius profile. It is seen that within the BL (z ^ (5„) the profiles match the theoretical 
solution very well, while in the region just outside the boundary layer where plume 
emissions occur, all measured profiles are generally less steep than the Prandtl-Blasius 
profile. This feature is also similar to that observed by (Zhou & Xia 20101) and will be 



L0| 
ft 



further discussed in Sec. |3.6[ On the other hand, it is seen from figure |5] that profiles 
obtained at different 9 seem to have different degrees of deviation from the Prandtl- 
Blasius profile. This can be seen more clearly in figure |6] where we show two examples 
in which profiles for different 9 but with values of Ra close to each other are plotted 
together along with the theoretical PB profile. This result suggests that the shape of the 
velocity profile near the plume-emission region is modified by the tilting angle. It is also 
noted that the profiles measured with 9 = 1.0° show strong deviations from the linear 
dependence with zero interception. We shall come back to this when discussing boundary 
layer scalings in the next section. 

Taking Umax as the characteristic velocity of LSC, we define the Reynolds number 
Re = UmaxH/v and plot Re as a function of Ra and for different 9 in figure [Tj^a). 
When fitting a power-law to the data for different 9 separately, they all produce an 
exponent close to 0.43. To better compare the amplitude of Re for different 6*, we fix 
the scaling exponent at 0.43 and fit power laws to the different data sets again. This 
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Figure 5. Profiles normalized by their respective maximum velocity Umax{Ro) and the corre- 
sponding viscous boundary layer thickness Sv{Ra) with tilt angles 6 = 0.5° (a), 1.0° (b), 2.0° 
(c) and 3.4° (d). The solid line in each plot represents the theoretical Prandtl-Blasius profile. 



gives Re = (0.185 ± 0.002, 0.182 ± 0.003, 0.206 ± 0.001, 0.203 ± 0.001) x Ra°-^^, where the 
amplitudes in the brackets are for 9 = 0.5°, 1.0°, 2.0°, and 3.4°, respectively. These results 
show that in general the values of Re with larger 9 are larger than those with smaller 
9. In an earlier study of the effect of cell titling, Ahlers et al. (2006) have found that 
Re obtained indirectly from temperature measurement increases with the tilted angle, 
which is consistent with the trend observed here. We note also that the value of the scaling 
exponent of Re obtained from many previous studies, and sometimes under nominally 
similar conditions, varies over a rather wide range from 0.43 to 0.55 (see for example. 
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for such variations is not completely clear at present. A detailed study on this issue is 
beyond the scope of this paper. For interested readers, we refer to Sun fc Xia| ( 2005 ) who 
offered an explanation that can account some of these dispersions in the exponent. 

From the measured profile of the RMS velocity (see figure [8|, we can define another 
Reynolds number Re^r = a„iaxH/v, which is shown in figure l7[b) as a function of Ra 
in a log-log scale for the four tilt angles. Here it is seen that Re^ does not seem to have 
an obvious dependence on 9. We therefore fitted a single power law to all four data sets 
on the plot, which gave Rca = O.OOT/Ja"-^^^""^. The value of the exponent is somewhat 



larger than 0.5 that was obtained from several previous studies (Xin et al. 1996 Xin 



& Xia 1997 Qiu & Xia 1998a|6 Sun et al. 2008). But given the uncertainties in the 
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Figure 6. Normalized profiles measured at different tilt angles 9 but with approximately the 
same value of Ra. Profiles in (a) have a nominal value of Ra = 5 x 10* and in (b) have a nominal 
value of Ra — 1.5 x 10®. In both figures the symbols are: inverted triangles (9 = 0.5°); squares 
{6 = 1.0°); triangles {9 = 2.0°); and circles {9 = 3.4°). The solid line in each plot represents the 
theoretical Prandtl-Blasius profile. 



experimental measurements, it is hard for one to attach too much significance to this 
difference. 

3.3. The viscous boundary layer and its scaling with Ra and Re 

We define the thickness Sy of the viscous boundary layer through the "slope-method" 
as shown in figure [s] where a mean velocity U{z) (circles) profile and the corresponding 
standard deviation profile cr^iz) (crosses) are shown, which are measured at Ra — 4.2 x 
10® with 9 = 3.4°. It is seen that (5„ is defined as the distance at which the extrapolation 
of the linear part of U{z) equals its maximum value Umax^ i-e- <5„ = Um.ax[dU / dz\z=Q\~'^ . 
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Figure 7. (a) Re based on the maximum horizontal velocity Umax, and (b) Re^ based on 
the maximum velocity fluctuation (Jmax as a function of Ra for different tilting angles. Inverted 
triangles: 6 = 0.5°; squares: 1°; triangles: 2°; circles: 3.4°. The dashed lines in (a) represent 
power-law fits to the individual data sets all with a scaling exponent —0.43 (see text for the 
fitting results). The solid line in (b) is a power law fit to all data sets in the plot, which gives 
Re. = 0.007J?a°-55±°«^ 




V ^ z (mm) 

Figure 8. Determination of the viscous boundary layer thickness (5„ through the slope-method 
from the mean horizontal velocity profile U{z) (circles) and the thickness 5„ from the standard 
deviation profile au (crosses) . The measurement was made near the bottom plate with tilt angle 
e = 3.4° and at Ra = 4.2 x 10®. 



A length scale 6a can also be defined from the profile of ct„(z) where cr„ reaches its 
maximum value. For the present example, the values for the two boundary layer length 
scales Sj; and 5a are found to be 4.20 and 6.05 mm, respectively. For ease of reference, 
the values of are listed in Table [T] 

We now examine the scalings of the boundary layer thickness with both the Rayleigh 
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Figure 9. Measured viscous boundary layer thickness 5^ (a) versus Ra and (b) versus Re for the 
four tilt angles: 9 — 0.5° (inverted triangles); 1.0° (squares); 2.0° (triangles); and 3.4° (circles). 
The dashed lines are power-law fits Sy/H = AiRa^^ and Sv/H = A2Re^'^ to the respective data 
sets, with the fitting results listed in Table [2] 



number Ra and the Reynolds number Re. In figures |9|^a) and (b) we plot the measured 
viscous boundary layer thickness 5^ vs, respectively, Ra and Re for the four tilt angles. 
The lines in the figures represent the best power-law fits Sv/H — AiRa^^ and 5v/H = 
A2Re^^ to the respective data sets and the obtained fitting parameters are listed in 
Table[2] Also shown in the Table for comparison are results obtained in cells with different 
geometries and using different methods. It is seen from the table that for small tilt angles 
{6 = 0.5 and 1.0°), the exponents are essentially the same and within the experimental 
uncertainties the i?e-scaling exponent may be taken as the same as that predicted for a 
Prandtl-Blasius boundary layer, i.e. Sy ~ i?e~^/^. For larger titling angles, there appears 
to be a trend for both (5i and /32 to decrease (absolute value increases) with increasing 
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Figure 10. The scaling exponent /32 versus cell tilt angle 8, where /32 is obtained from the 
power law fit ~ i?e^^. The dashed line indicates /32 = —0.5 for a Pradtl-Blasius laminar 

boundary layer. 



9. It thus appears that titling the cell by over 1° is a rather strong perturbation to 
the BL, at least for its scaling. The situation for the amplitude of viscous boundary 
layer thickness 5^, is a bit more complicated. From both figures |9];a) and (b) it seems 
that at lower values of Ra {Re) the BL thickness increases with increasing tilting angle, 
except for 6 = \° . For this latter titling angle, 5y appears to have an overall upward 
shift from the rest data sets. While we do not know the exact reason(s) for this, we note 
from figure [6] that the profiles for this tilt angle seem to have a nonzero intercept on the 
horizontal axis. This appears to suggest that the origin of the z-axis for this was somehow 
shifted. But even if this is the case, the relatively small "shift" cannot account for the 
large "deviation" of this 5y from the rest data sets (assuming there is indeed something 
"wrong" with this data set). Aside from the amplitude, the behavior of the Ra- and 
i?e-scaling exponents may be summarized as follows. For small tilting angle {9 ^ 1°), 
the effect of tilting is to lock the azimuthal plane of the LSC (or restrict its azimuthal 
meandering range) but the BL is otherwise not strongly perturbed and scaling wise the 
BL is approximately Prandtl-Blasius type. For relatively large titling angle {9 > 1°), the 
BL appears to be strongly perturbed as far as scaling is concerned and the magnitude 
of the scaling exponent increases with titling angle, i.e. the BL thickness Sy decays with 



increasing Ra (Re) with a steeper slope. The situation is illustrated in figure 10 where 
^2 is plotted as a function of the tilt angle 9. 

In addition to the boundary layer thickness Sy determined from the mean horizontal 
velocity profile, another length scale can also be defined based on the profile of the 
horizontal r.m.s. velocity cr„, which may be called the r.m.s. velocity boundary layer 
thickness, as defined in figure [8| In fi gures [TT|a) and (b) we plot 6cr versus Ra and Re^ 
respectively. The i?a-scaling exponent varies from —0.15 to —0.29, which appears to 
follow similar trend as that of i.e. its absolute value increases with increasing 9. But 
it and that of i?ecr-scaling exponent show significant difference with those obtained in 
previous studies. Table [2] shows the fitting results of 6c = A^Ra^^ and 6 a = A4_Re^*. 

Now we compare our result with previous experimental results obtained in the cells 
with different geometries. As shown in Table[2] the value of /3i obtained in both cylindrical 
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Table 2. Fitting results for the normalized viscous boundary layer thickness Sv/H determined 
from the mean horizontal velocity profile and S^/H determined from the rms horizontal velocity 
profile. The fitting parameters Ai and /3i (i = 1,2,3,4) are defined through the power laws: 
= ^liJa"!, S^/H = A2Re'^\ S^/H = A^Ra^'' , and S^/H = A^Re^^. The control param- 
eters Ra and Pr and cell geometry of measurements are also listed. Also show n i n the table 
are res ults from some previous experiment s. The sour ces are: a. |Xin et al. (1996 1; b.|La m et al.\ 

Jonnnl. _ \r\:.. c. v:„l /TTrino i /l.^4-4-^^\. J \r\:., c. v:„ / i nno !, \ /„: J II T IL'-.L _/ /onnoi. 



( 2002| ); c. |Qiu fc Xi"a| ([T998a ) (bottom); d. |Qiu fc Xi a (1998b) (side wall); e. |Sun et a l 
and f. preisent work. (Note: the cell tilt angle 6 is indicated as when it was not ment 



(2008), 

and f. present work. (Note: the cell tilt angle t) is indicated as when it was not mentioned in 
the respective papers and we assume the cell was nominally leveled in those cases.) 



10» 




a 

a. 



Figure 1 1 . Scalings of the boundary layer scale determined from the measured rms velocity 
profiles: (a) versus Ra and (b) versus Rcc- The symbols represent: 8 = 0.5° (inverted triangles), 
1.0° (squares), 2.0° (triangles), and 3.4° (circles). The lines are power law fits Sa/H = A^Ra'^^ 
and 5a/H ~ A^Re^^ to the respective data sets, with the fitting results listed in Table|2] 
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and cubic geometries and measured near the bottom plate of the cell is —0.16. In all these 
previous measurements, the Reynolds number based on the maximum horizontal velocity 
near the plate was also obtained and they gave a scaling exponent 7 = 0.5 via Re ~ i?a^. 
From this we obtain /32 = —0.32. In these studies, the convection cell was nominally 
leveled, i.e. not intensionally tilted. In the present study, for the small tilting angle cases, 
where we assume the BL is not strongly perturbed, the measured /3i ~ —0.19 when 
combined with combined with 7 = 0.43 give a /32 — —0.45 ± 0.04 (note that the actual 
value of /32 are obtained from fitting the Sy vs Re data, not from the relationship between 
the exponents). If we take these values to be close to the Prandtl-Blasius result, then 
scaling wise the viscous BL in a cylindrical geometry is also of a Prandtl-Blasius type, 



as was already found in a rectangular cell (Sun et al. 2008). For the relatively large 
deviations found in the untilted case, it may be attributed to the random azimuthal 
motion of the LSC. 

Finally we remark that as far as the scaling of the viscous BL is concerned, there is 
no theoretical prediction for the dependence of 6^ on Ra, only that on Re (for example, 
^i, ~ Re~^/'^ for the Prandtl-Blasius BL). In the literature, it is sometimes stated that 5^ 
should scale as Ra~'^/^ for the Prandtl-Blasius BL. This is based on the assumption that 
Re ^ Ra}/^. From above we have seen that the scaling exponent of Re with Ra varies 
over a rather wide range. It is therefore more meaningful to talk about the scaling of 5„ 



with i?e, rather than with Ra. We further note that in Sun et al. (2008j) it was found that 
5v ^ Ra^^'^'^ and Re ^ i?a°-^^, which together give Sy ^ Re~^'^^. In the present case, 
we have Sy ~ Ra^^-^ and Re ^ Ra° '^^, which together give Sy ^ Re~°-'^^'^° °^ . Whether 
this is fortuitous or there is something deep here remains remains to be explored. 

3.4. Fluctuations and statistical properties of the velocity field in the boundary layer 

In previous BL measurements in the cylindrical cell, owing to the nature of the dual- 
beam incoherent cross-correlation technique employed ( Xin et aL|1996 Lam et a/.|2002 ), 



only time-averaged velocity profiles are measured and no time-dependent quantities are 
obtained. It is therefore interesting to examine these quantities and compare them with 
similar quantities obtained in other type of turbulent flows. Figure [12] shows the time 
series of both the horizontal component u{t) (left panel) and the vertical component 
w(t) (right panel) of the velocity, measured at various positions from the plate. The 
corresponding velocity histograms are shown in figure |13| The measurements were made 
at Ra = 2.4 x 10* and 9 = 1°. We show the velocity trace at several typical positions: (i) 
inside the thermal boundary layer, (ii) around the thermal boundary layer, (iii) around 
the viscous boundary layer; (iv) at the position of the max:imum velocity; and (v) far 
away from the boundary layers. The figures show that the absolute horizontal velocity is 
much higher than vertical velocity at each position. One general feature we observed is 
that velocity time series and histograms look similar for different tilting angles. For this 
reason, we show here results for only one tilting angle. 

At Ra = 2.4 X 10*, the viscous BL thickness is 6y = 5.07 mm. It is seen from fig- 
ures 12 'a) and 12 'b) that at positions inside the BL, the horizontal velocity u{t) skews 
toward the positive side, i.e. the velocity is skewed toward the mean flow direction. This 
may be understood by the fact that close to the viscous sublayer the flow speed is very 
close to zero and a fluctuation smaller than the mean would mean a flow reversal, which 
is a rather rare event. Once outside of the BL, one observes more symmetric fluctuations 
around the mean velocity. For the vertical velocity w(t), its mean velocity is very small 
at most positions. But the fluctuation increases significantly when the position is outside 
of the BL, which are signatures of plume emissions at these positions. These properties 
can also be seen from the velocity histograms shown in figure [13) A notable difference of 
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Figure 12. Time traces of horizontal u{t) (left panels) and vertical w{t) (right panels) velocity 
components measured at Ra = 2.4 x 10® and 6 = 1°, at x = and different distances z from 
the bottom plate. 



the present results from those observed in Sun et at. ( 2008 1 is that for positions outside 



of the BL the horizontal velocity fluctuates more or less symmetrically around the mean, 
rather than skewed toward the negative as seen in the rectangular cell. 

The statistical properties of the velocity may be characterized more quantitatively 
by its root-mean-square (r.m.s.) value and its skewness, which are shown in figure 14 
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Figure 13. Histograms of (a) the horizontal velocity and (b) vertical velocity wit') 
measured at various distances from the bottom plate with Ra = 2.4 x 10* and 9 — 1°. 
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Figure 14. Profiles of (a) the normalized rms velocity {(7^) and of (b) the skewness Su {Sw) 
measured at Ra = 2.4 x 10* and 9 — 1°. The vertical distance z is normalized by the velocity 
boundary layer thickness 5^. In both plots the circles represent those for the horizontal velocity 
component u and the crosses represent those for the vertical velocity component w. 



Figure |14[ a) plots the velocity r.m.s cr„ and normalized by maximum horizontal 
velocity Umax versus the normalized distance z/S^v Figure |14| ^b) shows the skewness 
profiles Su = {{u - - {u))^)f/^ and Su, = {{w - - (u;))'))^/^ for 



the horizontal and vertical velocities, respectively. Similar to |Sun et al. (20081, our result 
could not tell whether cTu, favors a power law or a logarithmic scaling with the distance 
z, even though our measurement had a much higher spatial resolution. This is partly due 
to the limited size of the measurement area. 

3.5. Properties of shear stresses and near-wall quantities 

One of the advantages of PIV measurement is that it enables one to measure the hori- 
zontal and vertical velocities at the same time, so that one can calculate the Reynolds 
shear stress tji = —p{z){u'w'). Here u' and w' are the fluctuations of the horizontal and 
vertical velocity components respectively, p{z) is the z-dependent fluid density. Viscous 
shear stress is defined as Ty = fi{z)du/dz, where ij,{z) is the dynamic viscosity dependent 
on position z. The Reynolds stress represents the transport of momentum by turbulent 
fluctuations, whereas the viscous stress describes the momentum transfer by viscosity. 
The total shear stress is then r = t„ + r^j. 

Figure [TS] plots the profiles of the viscous shear stress, Reynolds stress and total stress 
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Figure 15. Viscous stress Tv (solid squares), the Reynolds stress tr (open circles) and the 
total stress r (open squares) as functions of the normalized distance from the plate for (a) 
Ra = 4.2 X 10® and (b) Ra = 9.9 x 10^ both with 9 = 3.4°. 
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Figure 16. The 7?a-dependence of (a) the wall shear stress r„, (6) the friction velocity Ur, (c) 
the wall thickness (viscous sublayer) 5^, and (d) the friction coefficient c/, with power law fits 
shown as solid lines. The symbols represent data for different tilting angles: 6 = 0.5° (inverted 
triangles), 1° (squares), 2° (triangles), and 3.4° (circles). 



for (a) Ra = 4.2 x 10*, and (b) 9.9 x 10*. It is seen that both Ra have the same qualitative 
features. Here the examples are for /3 = 3.4°, and results for other tilting angles are 
similar. Near the plate, it is seen that the Reynolds stress r^j is close to zero, while 
the viscous shear stress Ty is maximum because of the large velocity gradient du/dz 
at the wall. So the total stress at the wall r^(= t(0)) comes almost entirely from the 
contribution of the viscous shear stress. Moving away from the plate, the velocity gradient 
becomes smaller and the viscous shear stress decreases to zero. The Reynolds stress tr 
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Figure 17. The TJe-dependence of (a) the wall shear stress r^, (b) the friction velocity Ut, (c) 
the wall thickness (viscous sublayer) 5w, and (d) the friction coefficient c/. The symbols represent 
data for different tilting angles: 9 = 0.5° (inverted triangles), 1° (squares), 2° (triangles), and 
3.4° (circles). Power law fits are indicated in the figure. 
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Figure 18. Measured horizontal velocity profiles normalized by wall units for four different 

Ra all taken at tilting angle 9 — 3.4°. 
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increases and attains its maximum at z « \.^5y. It then decreases to around zero z « 25^ 
and becomes negative in the bulk flow. It is also seen clearly from the figure that r/j 
and Ttj cross at z l.bSy, where tji is close to its maximum value. This suggests that 
the momentum transfer in the outer region is dominated by turbulent fluctuations. But 
in the viscous boundary layer, the momentum transfer is still dominated by the viscous 
diffusion, which implies that the viscous boundary layer is still laminar in this range of 
Ra. 

With the measured near-wall high-resolution velocity field, we are now in a position 
to check the dynamic wall properties in turbulent thermal convection. We first consider 
the scaling of four basic wall quantities with both Ra and Re. These are the wall shear 
stress Tw, the skin-friction velocity Ur = (r/po)^^^, the viscous sublayer length scale 
5w — vq/ut, and the skin-friction coefficient c/ = t/ poU^^^. Here po = p{z — 0) and 
uq = v{z = 0). Figure 16 shows the scaling of these quantities with Ra. It is seen that 
within experimental uncertainties there is no difference between data with different 9. 
This suggests that tilting the cell does not have any appreciable effect on BL properties 
near the wall. Without differentiating the different data sets, power law fits to all data 



yield t„, ~ Rd 



0.63 



Ra: 



0.32 



Ra 



-0.37 



and Cf ^ Ra 



-0.19 



In a rectangular cell. 



Sun et al (2008) found for the same quantities the fitted power law exponents 0.86, 0.44, 
—0.50, and —0.28 respectively. It is seen that the absolute values of these exponents 
are all larger than those obtained in the present experiment. There is no theoretical 
prediction for the iZa-scaling of these quantities in turbulent thermal convection, so we 
do not know what the difference means. 

It will be more useful perhaps to examine the scaling of these quantities with the 
Reynolds number i?e, since theoretical predictions exist for such scalings for wall-bounded 
shear flows ( Schlichting fc Gersten|[2000 ) . Figure 17 plots these quantities as a function 
of Re, the symbols are the same as in figure 16 For the quantities t^, "r, ^w, and Cf our 
results give the exponents 1.46, 0.75, —0.86, and —0.46. For a laminar boundary layer over 
a flat plate, the theoretically predicted 'classical' exponents for these quantities are 3/2, 
3/4, —1, and —1/2 respectively. One sees that within the experimental uncertainties there 
is an excellent agreement between the present experiment and the theoretical predictions 
for all the wall quantities except for (5^,, which is a bit smaller. For reference, the previous 
measurement in rectangular cell gives 1.55, 0.8, —0.91, and —0.34 for the corresponding 
quantities ( Sun et a7^|2008 ). 

To further compare the present system with classical boundary layers, we examine ve- 
locity profiles in terms of the wall units. Figure [18] shows the normalized mean horizontal 
velocity profiles for four different values of Ra taken at 6* = 3.4° in a semi-log plot, here 
u"*" = u[z)/ut and z+ = z/5w. The linear scaling of over z^ in the viscous sublayer 
below z+ < 5 is refiected quite well by the measured profiles confirming that the bound- 
ary layer is not turbulent in the present range of Ra and Pr. The velocity normalized 
by wall unit decrease after reaching the maximum value in z"*" ~ 10. Comparing to the 
same quantity measured in the rectangular cell (Sun et al. 2008), however, our result 



shows some deviation from the theoretical profile. This is a reflection of the fact that in 
the cylindrical cell it is more difficult to measure the profile accurately very close to the 
wall. 



3.6. Dynamical scaling and the shape of velocity profiles in the boundary layer 



The dynamic scaling method of Zhou & Xia (2010) has been found to work well when 



tested in quasi-2D experiment and 2D numerical simulations (Zhou et al. 2010 2011) 



But it has not been examined in 3D experiments. Here we investigate how it works in our 



cylindrical geometry. As the method has been well documented elsewhere (Zhou & Xia 
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Figure 19. Comparison between profiles obtained in the dynamical frame (u*{z*), red circles) 
and the laboratory frame {u[z), blue squares), measured at different tilt angles 6 but with 
comparable values of Ra. Also shown for comparison is the theoretical Prandtl-Blasius laminar 
velocity profile (solid line), (a) 6 = 0.5°, Ra = 5.77 x 10^ (b) 6 = 1.0°, Ra = 6.00 x 10*; (c) 
e = 2.0°, Ra = 5.85 x 10*; and (d) 6 = 3.4°, Ra = 6.69 x 10*. 



2010 Zhou et al. 2010 2011), we will only give a brief description of it here. From the 
measured instantaneous velocity profile u{z,t) one can obtain an instantaneous viscous 
boundary layer thickness Sjj{t) using the same 'slope' method as used for the mean 
velocity profiles. A local dynamical BL frame can then be constructed by defining the 
time-dependent rescaled distance z* (t) from the plate as 



z*{t) = z/S,it). 



(3.1) 



The dynamically time averaged mean velocity profile u*{z*) in the dynamical BL frame 
is then obtained by averaging over all values of u(z,t) that were measured at different 
discrete times t but at the same relative position z* , i.e. 



i*{z*) = {u{z,t)\z — Z*Sv{t)). 



(3.2) 



Figure 19 shows the mean velocity profiles measured in the laboratory and the dy- 
namical frames respectively at the four tilting angles and for comparable values of Ra (as 
indicated in the figure caption). These results show that the dynamical scaling method 
appears to be more effective for larger values of 9. This may be understood based on the 
fact that a larger tilt angle will place stronger restriction on the azimuthal meandering 
of the LSC so that it has less fluctuations in the horizontal direction perpendicular to 
the mean flow. We note, however, regardless of the tilt angle, the method works less 
effectively than it is in quasi-2D experiment and 2D simulations. 
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Figure 20. The shape factor H — 5d/Sm of profiles u*{z*) (red circles) obtained in the dy- 
namical frame and of u{z) (blue squares) obtained in the laboratory frame as a function of Ra 
and for different titling angles, (a) 9 = 0.5°, (b) 1.0°, (c) 2.0°, and (d) 3.4°. The dashed line 
represents the value of 2.59 for the theoretical Prandtl-Blasius laminar BL. 



A more quantitative approach to characterize the shape of the mean velocity profiles is 
to investigate their shape factor H — Sy/Sm defined as the ratio between the displacement 
thickness Sy and the momentum thickness Sm, where 



u{z) 



]dz, and 6„ 



u{z) u{z) 



dz. 



(3.3) 



Since u(z) decays after reaching its maximum value, the above integrations are evaluated 
only over the range from 2; = to where u(z) = Umax - For our profiles the obtained shape 
factors range between 1.9 to 2.3, which are smaller than H = 2.59, the value for a laminar 
Prandtl-Blasius boundary layer. A shape factor smaller than the theoretical value means 
the corresponding profile will approach its asymptotic value (the maximum velocity) 
slower than the theoretical profile does. 

In figure [20] we show the shape factor H for mean velocity profiles obtained in the 
laboratory and dynamical frames respectively for the four tilting angles and for all Ra 
measured. The dashed lines in the figure indicate the Prandtl-Blasius value of 2.59. 
It is seen that, despite the data scatter, there is a general trend that for both lab- and 
dynamical-frame profiles the deviation from the Prandtl-Blasius profile increases with Ra. 
This is no surprise, since, as the convective fiow above the BL becomes more turbulent 
with increasing Ra, the BL itself will experience stronger fluctuations and hence larger 
deviations from the laminar case. This finding that the dynamical rescaling method works 
better for smaller Ra than larger ones is consistent with those found in DNS studies in 



the same geometry by Stevens et al. (2012) for the temperature profile and by Shi et al. 
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Figure 21. Examples of instantaneous horizontal velocity profiles scaled by the instantaneous 
kinematic BL thickness and the instantaneous maximum velocity (measured at Ra = 6.00 x 10* 
and 9 = 1°). The corresponding instantaneous shape factor is also indicated on the plot. The 
solid curves are the Prandtl-Blasius velocity profiles. 



(2012) and by Scheel et al. (2012) for the velocity profile. The second feature is that for 
all 9 and Ra the profiles obtained in the dynamical frame in general show some degree 
of improvement towards that of Prandtl-Blasius value as compared to those obtained in 
the laboratory frame. We also note that the "degree of improvement does not seem to 
have an obvious dependence on Ra, which is also consistent with the findings of [Zhou 



Xia (2010); Zhou et al. (2010). 



Some insight can be obtained by examining the rescaled instantaneous velocity profiles. 
Figure [21] show examples of rescaled instantaneous velocity profiles, where the distance 
from the plate has been normalized by the instantaneous BL thickness corresponding to 
that moment and the velocity has been normalized by the instantaneous maximum hori- 
zontal velocity. It is seen that there are quite few cases where the rescaled instantaneous 
velocity profile is rather close to the theoretical Prandtl-Blasius profile (up to the point 
of the maximum velocity) and deviations of the instantaneous shape are likely caused 
by distubances such as plume emissions. Also shown in the figure are the shape factor 
H{t) of these instantaneous profiles. To quantify how the instantaneous profiles are dis- 
tributed with respect to the Prandtl-Blasius profile, we examine the PDF of the shape 

where H^^ 
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plots the PDFs of 



factor difference SH{t) = H{t) - H^^ where H^" = 2.59. Figure 
5H{t) for the 4 tilting angles and for all measured Ra respectively. Despite the seem- 
ingly large variations among them, these PDFs show the general trend that the rescaled 
instantaneous profiles measured at lower values of Ra (;$ 1 x 10^) are more of the time 
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Figure 22. PDFs of the shape-factor difference 5H between those of the rescaled instantaneous 
profiles and that of the Prandtl-Blasius profile, measured at (a) 6 = 0.5°, (b) f .0°, (c) 2.0°, and 
(d) 3.4°. 



having a shape closer to that of the Prandtl-Blasius profile and that for higher values of 
Ra the peak of the PDFs shift to smaller values of H . This indicates that with increasing 
Ra the profiles around the BL thickness becomes more rounded, i.e. the approach to the 
maximum velocity becomes slower and slower. We further note that these general trends 
are true across all tilt angles. Another feature observed in the present 3D case is that 
we did not find any strong correlation between the instantaneous BL thickness 5^ and 
the velocity u{t) just above the BL. This is in contrast to the finding in the quasi-2D 
experiment where 5„ and u(t) are found to have a strong negative correlation, i.e. a large 
velocity above would exert a stronger shear and therefore thins the BL thickness (jZhou] 
I&: Xia| 2010 1. This result suggest that in certain aspect the BLs in the 3D and in the 
2D/quais-2D cases are dynamically different. 



4. Summary and conclusions 

We have conducted an experimental study of velocity boundary layer properties in 
turbulent thermal convection. High-resolution two-dimensional velocity field was mea- 
sured using the particle image velocimetry (PIV) technique in a cylindrical cell of height 
H = 18.6 cm and aspect ratio close to unity, with the Rayleigh number Ra varying from 
10® to 6 X 10^ and the Prandtl number Pr fixed at ~ 5.4, with the convection cell tilted 
with respect to gravity at angles = 0.5°, 1°, 2°, and 3.4°, respectively. Measurements 
made with small 9 are aimed at studying BL properties under more steady shear, but 
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the BL itself is assumed to be unperturbed otherwise. For large values of 9 we wish 
to examine how the BL responds to relatively large perturbations. We also examined 
effectiveness of the dynamical BL scaling method in a three-dimensional system. 

It is found that the Reynolds number Re (= UmaxH/v) based on the maximum mean 
horizontal velocity scales with Ra as Re ~ Ra^-^^ and the Reynolds number Re„ (= 
c^maxH/v) based on the maximum rms velocity scales with Ra as Re„ ^ Ra'^'^^. Both 
exponents do not seem to have an apparent dependence on the tilt angle. On the other 
hand, the amplitude of Re seem to show a weak increasing trend with 6. 

With the measured horizontal velocity, we obtain two length scales, i.e. the viscous BL 
thickness Sy based on the mean horizontal velocity profile and the length scale S„ based 
on the rms horizontal velocity profile. It is found that as far as scaling with the Reynolds 
number Re is concerned, the behavior of Sy can be divided into two regimes according to 
the tilting angle of the cell. For 9 < 1°, it is found that dy ~ Re-°-'^'^'^°-°^ , which within 
experimental uncertainty may be considered to be consistent with that of the Prandtl- 
Blasius BL. It thus appears that the main eifect of tilting the cell is to restrict the 
azimuthal meandering of the large-scale circulation but the BL is otherwise not strongly 
perturbed. For 9 ^ 1°, the absolute value of the exponent is found to increase with 9 
and in this case the BL may be considered to be strongly perturbed. It is found that 
the scaling exponent of Jo- with respect to Ra {Re) does not have a strong dependence 
on 9 as 6v does. But similar to Sy, the absolute values of these exponents increase with 
increasing 9. 

It is also found that tilting the cell modifies the velocity profile in the BL region, i.e. 
for different tilt angles the shape of profiles is diff'erent. But for the same tilting angle the 
velocity profiles measured at different Ra can be brought to collapse on a single curve 
when the mean velocity is normalized by the maximum velocity Umax and the distance 
from the plate by the viscous BL thickness 6^. 

With simultaneously measured horizontal and vertical velocity components, we also 
obtain the Reynolds stress tr in the velocity boundary layer. It is found that tr is stronger 
in the mixing zone comparing with the rectangular cell. The wall quantities such as the 
wall shear stressr^, the viscous sublayerJ^, the friction velocity Ur are also measured. 
Their scaling exponents with the Reynolds number are very close to those predicted 
for classical laminar boundary layers, which is also consistent with the measurement in 
rectangular cell. 

Regarding the dynamical scaling method, we found that the method in general works 
better when the cell is tilted at larger angle 9 than it does at smaller angles, but the 
effect is somewhat marginal. With respect to the influence of Ra, it is found that in 
general the mean velocity profile sampled in both the laboratory and dynamical frames 
are more closer to the Prandtl-Blasius profile at smaller values of Ra than they are at 
larger Ra, which is consistent with findings from previous DNS studies. Moreover, it is 
found that for smaller values of Ra (< 1 x 10^) the PDF's of the shape factor H for 
the rescaled instantaneous profiles exhibit a peak close to that for the Prandtl-Blasius 
profile, whereas for larger vahies of Ra the peaks shift to smaller values of H , indicating 
the profile's approach to the maximum velocity becomes slower and slower with increasing 
Ra. Another finding is that the effectiveness of the dynamical scaling method, in terms 
of its ability of bringing the mean velocity profile closer to that of Prandtl-Blasius profile, 
does not have any apparent dependence on Ra. Our general conclusion is that as far as 
the effectiveness of the dynamical scaling method is concerned the influence of titling 
angle is much smaller than that of the Raylcigh number Ra. We note that the Prandtl- 
Blasius boundary layer theory is a 2D model, so it is perhaps no surprise that the dynamic 
method works less well in 3D than in 2D. 
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